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The spin-1 Ising (BEG) model has been simulated using a cellular au- 
tomaton (CA) algorithm improved from the Creutz cellular automaton (CCA) 
for a face-centered cubic (fee) lattice. The ground state diagram {k, d) of 



Ch ' the fee BEG model has ferromagnetic (F), quadrupolar [Q) and staggered 
O ' 

, quadrupolar {SQ) ordering regions. The simulations have been made in 

\Q • the staggered quadrupolar region for the parameter values in the intervals 
> 

ON : -24 < = L)/J < and -3 < A; = K/J < . The phase diagrams 



on the {kTc/J, d) and the {kTc/J, k) planes have been obtained through 
k = —3 and d = —4 lines, respectively. The staggered quadrupolar order- 
ing region separates into five ordering regions (A^B^a), A^B^f), AB (type-1), 
AB{tjpe-lT) and AB^lf)) which have the different stoichiometric Cu-Au type 
structures. 

Key words: Spin-1 Ising (BEG) model; Cellular automaton; phase dia- 
gram; face-centered-cubic lattice 
1 . Introduction 

The Blume-Emery-Griffiths (BEG) model [1] is a spin-1 Ising model with 
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very rich and interesting phase structure. Some of models in the statistical 
mechanics can be regarded as special cases of a general Ising model. The 
lattice gas Hamiltonian may be written in terms of site occupation operators 
Pi" constituting by the Ising spin energy [2 — 12]. In this case, the generalized 
Hamiltonian is defined as [6, 13] 

^ = E E ^^P^P^ + E E + ^0 (1) 

where the I^J^ describes the interaction between species A and u on sites i 
and j, respectively. Second term in Eq. (1) describes the binding energy. 
The irrelevant terms are included in Hq. If site i is occupied by species A, 
Pi' is equal to 1, and zero otherwise. This energy definition is also valid for 
the ternary alloys [6]. This model can be related to spin-1 Ising model by 
the following transformations: 

Pf = SiiSi - l)/2 

Pf = S,{S, + l)/2 

The three components of the spin variables Si are related with each species 
of atoms: Si = -1 (A), +1 (B) and (C) for the ((AB)i_j;C2x) ternary alloys 
or Si = (Cu) and ±1 {Au) for Cu-Au type alloys. The Hamiltonian is then 
equivalent to the spin-1 Ising Hamiltonian 



Hi^-jJ2 S^S,-K S^S]+L ^ {SiS]+SfSj)+D ^ S^+h ^ Si+H^ 

<ij> <ij> <ij> i i 

(2) 
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where 



J = -{Jaa + Jbb)/4: - Jab/2 (3) 

K = -{Jaa + Jbb)/4: + Jab/2 - Jac - Jbc + Jcc (4) 

L = {Jaa - Jbb)/^ + {Jac - Jbc)/2 (5) 

D = z{2Jcc - Jbc - Jac) -Vb-Va + 2Vc (6) 

h = z{Jbc - Jac) ^Vb-Va (7) 

J , K , L, D , h and Hq are bilinear , biquadratic , dipol-quadrupol interaction 
terms, single-ion anisotropy constant, the field term and irrelevant terms, 
respectively. < ij > denotes summation over all nearest-neighbor {nn) pairs 
of sites. 

The spin-1 Ising (BEG ) model was firstly studied to determine the phase 
separation and superfiuid ordering in He^-He^ mixtures [1]. Variants of the 
model have been extensively studied because of the fundamental interest in 
the tricritical, reentrant and multicritical phenomena [14 — 16] of physical 
systems such as solid- liquid-gas systems [2], the multicomponent fiuids [3], 
the microemulsions [4], the metastable alloys [5,6, 17] and the binary alloys 
[18]. The BEG model has been investigated using many different approx- 
imate and simulation techniques in the special phase regions for difi^erent 
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lattice type and dimensions. These studies show that the fee BEG model 
with repulsive biquadratic coupling has the staggered quadrupolar (SQ) or- 
dering phase region in the {k = K/J, d = D/J) ground state diagram for 
the A; < — 1 parameter region [19 — 27], but there are the different results for 
the phase boundaries between the SQ ordering region and the other ordering 
regions. Furthermore, the SQ ordering region includes the Cu-Au {A^B, AB 
and AB^) type stoichiometric structures as the ground state ordering [24, 25]. 
The special orders of the Cu-Au {A^B, AB and AB^) type binary alloys are 
investigated by the MC [28 - 32], the CVM [33, 34] and the mean field the- 
ory [35,36] considering also the spin-1/2 Ising model, the effective medium 
theory (EMT) [31, 36] , the Bozzola Ferrante Smith (BFS) [32] method and 
the numerical calculation [37]. The Cu-Au type structures are also studied 
experimentally by the x-ray diffraction [39 — 42], the nuclear magnetic reso- 
nance (NMR) [43], the galvanic cell study [44], the electrical resistivity [45] 
and the electron diffraction [46, 47] in the binary alloys. Some of the theo- 
retical [31,36,37] and the experimental studies [39,41 — 47] for CuAu type 
structures have indicated that the AB ordering structure separated into two 
different structure as the modulated AB (type-II) phase and the AB (type- 
I) phase. The modulated AB (type-II) phase is produced from unit cells 
of the modulated AB (type-I) phase and is characterized by the existence 
of antiphase boundaries [46]. In this work, we have been investigated the 
Cu-Au type structures using an alternative algorithm improved from cellular 
automaton for the fee the BEG model with repulsive biquadratic coupling. 
In the previous papers, the Creutz cellular automaton (CCA) algorithm and 
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improved versions have been used successfully to study the properties of 
the critical behaviors of the Ising model Hamiltonians [15, 16, 48 — 67]. The 
CCA algorithm, which was first introduced by Creutz [53], is a microcanon- 
ical algorithm interpolating between the conventional Monte Carlo and the 
molecular dynamics techniques for Ising model. 

The basic aim of this study is to investigate the existence of the modu- 
lated AB (type-II) phase and the critical behavior of the Cu-Au {A^B, AB 
and ABs) type stoichiometric structures on the fee BEG model using the 
heating algorithm [15, 16] improved from cellular automaton. For this pur- 
pose, the fee BEG model is simulated on a cellular automaton (CA) through 
the k = —3 and d = —4 lines in the —24 < c? < and —3 < k <0 parameter 
region, respectively. The temperature variations of the four sublattice order 
parameters (ma, qa), the four sublattice susceptibilities (Xa), the lattice spe- 
cific heat {C/k) and the lattice Ising energy (Hj) were computed on the fee 
lattice with linear dimension L= 9 for J > 0, L — and H — 0. The finite 
lattice critical temperatures are estimated from the maxima of the specific 
heat (C/k) and the sublattice susceptibihties (Xa) to get the phase diagrams. 
The CA results for the 5% heating rate nearly agree with the CVM and MFA 
results for the critical lines. In addition, it is seen that the fee BEG model 
on the cellular automaton exhibits the modulated AB (type-II) phase as in- 
dicated by the theoretical studies [31, 36, 37] and the experimental results 
[39, 41 — 47] for Cu-Au type structures. 

2. Results and discussion 

The simulations have been performed on the face-centered cubic lattice 
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for linear dimension L— 9 with periodic boundary conditions using heating 
algorithm (The total number of sites is N = 4L^) . This algorithm is reahzed 
by increasing of 5% in the kinetic energy (Hk) of each site [15]. Therefore, 
the increasing value per site of Hk is obtained from the integer part of the 
O.OSiffe. In this study, +1 and —1 values of spin variable Si are related with 
species B (Au), while species A (Cu) is represented with value of S^. 
The computed values of the thermodynamic quantities are averages over the 
lattice and over the number of time steps (1.000.000) with discard of the first 
100.000 time steps during which the cellular automaton develops. 

The ground state diagram of the fee BEG model is illustrated in Fig. 1. 
The ground state phase diagram of the fee BEG model has the ferromagnetic 
(F), the staggered quadrupolar (SQ) and the quadrupolar (Q) ordering re- 
gions on the {k, d) plane. Furthermore the staggered quadrupolar ordering 
region separates into three regions which have the A^B^ ABand Ai?3orderings 
[24]. The initial simulations obtained for some dand /^parameter values in 
the interval —24 < d < Othrough the k — — 31ine are shown that the lattice 
order parameters (M, Q) do not have any sign in low temperature region 
while the behavior of the lattice Ising energy {Hi) indicate a phase transi- 
tion. Therefore, the sublattices have been followed to distinguish the type 
of the phases and the phase transitions. The temperature dependence of the 
sublattice order parameters are different each other at the ferrimagnetism 
and the antiquadrupolar phases [21,22,24]. 

The fee lattice can be built from the four interpenetrating simple cubic 
(sc) lattices, called sublattices. The sublattice order parameters (ma, ^q,) are 
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calculated as 



1 



(8) 
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a 
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where a indicate sublattices {a = a,b,c,d) (Fig. 2). 

According to the values of the sublattice order parameters, the model has 
the five different phases (P, F, A3B{a), AB3{f), AB (type-1)) at absolute 
zero temperature and two different phases {A^B^f ) and AB (typc-II)) above 
absolute zero temperature for —3 < A; < — 1 [21,24]: 

Paramagnetic (P) — {a — a,b, c, d), Qa — Qb — Qc — Qd > 0, 
Ferromagnetic (P) rUa ^ mi, ^ rric ^ ma 0, ^ q^, ^ q^^ > 0, 
Antiquadrupolar ^435(0) rria = {a = a,b,c,d), qa > Qb = Qc = Qd ^ 0, 
Ferrimagnetic A^Blf) rua > mb = rric = md> 0, qa > Qb = Qc = Qd > 0, 
{nia ^ qa, Oi^ a,b,c,d), 

Ferrimagnetic AB^lf) rUb ^ m^^ md> ma > 0, qb ^ qc ^ Qd >(la > 0, 
{iTLa = qa, a = a, 6, c, d), 

Ferrimagnetic AB (type-I) rUa = nib = I, mc = ma = 0, qa = qb = ^, 
qc^qd^ 0, {nia ^qa, a, b, c, d). 

Ferrimagnetic AB (type-II) rric — m^ > nia — mb > 0, qc — qd > qa — 
qb > 0, {nia = qa, = O) b, c, d). 

These phase definitions except AB ordering structure are compatible with 
the MFA and the CVM studies [21,24,25]. In this study ferrimagnetic AB 
ordering structure appears as AB (type-I) and AB (type-II) ordering struc- 
tures while the MFA and the CVM studies do not exhibit this separation. In 
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the AB (type-I) ordering structure, two of the sublattices are occupied by B 
{Si = +1 or —1) species ( = mi, = 1, qa = qb = 1), while the other sublat- 
tices are occupied by A {Si = 0) species (m^, = = 0, qc = qd = 0) on the 
fee lattice. In the ferrimagnetic AB (type-II) ordering structure, the A and 
B chains interchange periodically on the fee lattice. This is a modification 
of the AB (type-I) ordering structure [36, 46, 68 — 70]. 

2.1. The phase transitions on the Cu-Au type orderings 
The temperature variations of the sublattice order parameters {ma,qa), 
the sublattice susceptibility {xj, the lattice Ising energy {Hj), and the lat- 
tice specific heat {C/k) are calculated to determine the types of the phase 
transitions in the interval — 24 < d < through the A; = — 3 line and in 
the interval —S<k<0 through the d — —A line. For producing the phase 
diagrams, the finite lattice critical temperatures are estimated from the max- 
ima of the specific heat {C/k) and the sublattice susceptibilities (x^) on the 
lattice with L= 9 . 

In the interval — 8 < d < through the k — —?> line, the model has 
the paramagnetic A:iB{P) ground state ordering. In the A'iB{P) ordering 
structure, three of the sublattices are occupied by A {Si = 0) species ( 
rrib — rric — rud = 0, qh — qc — (ld — ^)-i while the other sublattice is occupied 
randomly by B {Si — +1 and —1) species {iria — 0, qa — 1)- The temper- 
ature dependences of the thermodynamic quantities are given in Fig. 3 for 
representing the A3B{a) — P phase transition at a selected {d = —4, k = —3) 
point on the {k, d) plane. It can be clearly seen that, as temperature in- 
creases, the sublattice order parameter qa decreases sharply from unity to a 
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value ( ~ 0.25) of the paramagnetic ordering and coincide with gr^, and 
qd at a transition temperature kTt/ J. Furthermore, the sublattice order pa- 
rameters (ma(a = a, 6, c, d)) do not exhibit any change from initial value. 
For this transition, the sublattice order parameters qa and the Ising energy 
are in s-shape nature (Fig. 3(a) and 3(c)) while the sublattice susceptibilities 
and the specific heat show a sharp peak at the transition temperature kTt/ J 
(Fig. 3(b) and 3(d)). The ^43^(0) — P phase transition is of first order for 
d — —4 and k — —3. The model has a similar behavior for some d parameter 
values in — 7 < d < parameter region through the A; = — 3 line. However, 
the phase transition is usually of second order for d parameter values in this 
region. On the other hand, the lattice order parameters (M, Q) indicate the 
paramagnetic ordering and they do not show any sign at the ^43^(0) — P 
phase transition in this parameter region (Fig. 3(a)). This result shows that 
the observing of the sublattices is necessary to decide the type of ordering 
and the phase transition. In the interval —7.3 < d < —7 for k = —3, the 
model exhibits the second order successive ^43^(0) — F— P phase transi- 
tion. As it is seen in Fig. 4, the sublattice order parameters continuously 
increase to a value of ferromagnetic F ordering (m^ = mb = rric = md = 0.20, 
Qb — Qc — Qd — 0.25) from the ^43^ (a) ground state ordering near the crit- 
ical temperature kTci/J. At high temperatures, the F ordering tends to 
paramagnetic P ordering. In this interval, the sublattice order parameters 
and the lattice Ising energy appear continuous near the critical temperatures 
kTci/ J and kTc2/ J while the lattice susceptibility an the sublattice suscep- 
tibility and the lattice specific heat have two characteristic peaks at these 
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temperatures. The all phase transitions are of the second order in this region. 

On the {k, d) plane, the model does not have any ground state ordering 
between the A^Bla) and AB (type-I) ground state ordering regions. How- 
ever, the CA results show that there are successive A^B{P)—A^B{f) — F — P 
phase transitions in the interval — 8 < d < — 7 through the k — —?> hne. The 
temperature variations of the thermodynamic quantities are illustrated in 
Fig. 5 for selected d = —7.7 parameter value. For this parameter, the 
model has AzB{a) ground state ordering . With increasing temperature, the 
A'iB{P) ordering tends to the ferrimagnetic A2,B{f) ordering near the critical 
temperature as CVM result [24]. At a certain temperature above the criti- 
cal temperature kTci/ J, the sublattice order parameters coincide {nia = 0.3 
, Qa — 0.35 {a — a, b, c, d)) with each other. In this case, the amount 
of the A (Si — 0) and the B (Si — +1 or —1) species on the each sub- 
lattices are equal and the ordering of the system changes from A^Blf) to 
F at the critical temperature kTc2/J- At a critical temperature kTcs/J 
above the critical temperature kTc2/ J, the model exhibits the second order 
F-P phase transition. As it is seen in Fig. 5, the sublattice susceptibihties 
and the lattice specific heat have the three characteristic peaks while the 
sublattice order parameters and the lattice Ising energy exhibit continuity 
for the successive phase transitions. This case indicates that the successive 
A^B{a)—A2,B{f) — F — P phase transitions are of the second order. 

For the interval — 16 < c? < —8 through the k = —3 line, the model 
takes the AB (type-I) ordering structure as ground state, thus — nib — 
1) Qa — Qb — 1, TTT'c — 'f^d — and Q'c = = at kT/J — 0. In Fig. 
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6, temperature variations of the thermodynamic quantities are illustrated 
for the d = —12 value through the A; = — 3 hne. It can be clearly seen 
that as temperature increases nia = mi, and Qa = Qb decreases sharply from 
unity to a value of the ferromagnetic ordering and coincide with nic — ma 
and qc = at a certain temperature, kTt/J. There are the successive AB 
(type-I)— F — P phase transitions. At the AB (typc-I)— F phase transition, 
the sublattice order parameters and the lattice Ising energy are in s-shape 
nature at the transition temperature kTt/J (Fig. 6(a) and 6(c)) and the 
sublattice susceptibility and the lattice specific heat have a sharp peak at 
kTt/J. For the F — P phase transition, the sublattice order parameters and 
the lattice Ising energy appear continuously near the critical temperature 
kTc/J (Fig. 6(a) and 6(c)). At the same time, the sublattice susceptibility 
and the lattice specific heat have a characteristic peak at kTc/J for this 
transition (Fig. 6(b) and 6(d)). The results show that the AB (type-1)— F 
and the F — P phase transitions for c? = — 12 parameter value are of the first 
order and of the second order, respectively. However, the 74S(type-I)— F 
phase transition is of the second order for the all values of the d parameter 
except d = —12 in the interval —15.5 < d < —8. 

In the intervals -15.8 <d< -14, -14 < ci < -12, -12 < ci < -10 and 
—10 < d < —8.3 of the phase space, the model has the ferrimagnetic AB 
(type-II) ordering between ferrimagnetic AB (type-I) and F orderings. The 
AB (typc-II) ordering is determined from the behavior of the sublattice or- 
der parameters ( m^, Qa) and peaks of the sublattice susceptibilities {Xa) and 
the lattice specific heat (C/k). For a selected d — —14.5 value, the temper- 
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ature variations of the thermodynamic quantities are shown in Fig. 7. The 
temperature dependences of the thermodynamic quantities show that the 
model exhibits the successive AB (type-I)— F — AB (type-II)— F — P phase 
transitions in these intervals. At the low temperature region of the succes- 
sive phase transition, the sublattice order parameters nia — irib — Qa — Qb 
decreases from unity to a value (~ 0.58) of the ferromagnetic ordering and 
coincide with rUc = rud = Qc = Qd parameters at the kTci/J temperature. 
This case indicates the phase transition from AB (type-I) to F order. At 
a critical temperature kTc2/J above the kTci/J, the sublattice order pa- 
rameters change their positions and part equally from = ~ 0.58 to 
produce the AB (type-II) ordering. It can be seen from Fig. 7(a) that the 
model exhibits AB (type-II) ordering in the certain temperature range. As 
the temperature increases, the sublattice order parameters coincide again 
with each other and there occur the F ordering at the critical temperature 
kTca/ J above the critical temperature kTc2/ J- Finally, the sublattice order 
parameters rria decrease to the value {rria ~ 0) of P ordering near a high 
temperature kTc/^j J above kTczj J ■ The sublattice order parameters (ma, 
qa) and the lattice Ising energy [Hi) show continuity (Fig. 7(a) and 7(c)) 
at the these successive phase transitions while the sublattice susceptibilities 
(Xq,) exhibit the broad peak with a strong peak at kTcxj J and a shoulder at 
kTc2/ J-, a weak peak at kTczj J and a characteristic peak at kTc^j J for the 
{).0'oHk value of the heating rate (Fig. 7(b)). However, the lattice specific 
heat {C/k) have a characteristic peak {kTci/ J ) , a broad peak [kTca/J 
and kTc2/J) and a strong peak at high temperature region {kTci/ J) (Fig. 
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7(d)). Therefore, all of the successive AS {type-T)-F - AB (type-II) -F - P 
phase transitions are of second order in these parameter regions. The simu- 
lations arc repeated at the different heating rates and the linear dimensions 
for checking the existence of the AB (type-I) —F — AB (type-II) phase tran- 
sitions at the different d values. The results of the sublattice susceptibihties 
and the lattice specific heat for the O.OIHk and 0.05Hk values of the heat- 
ing rate are illustrated in Fig. 7(b) and (d). It can be clearly seen from 
figure that the lattice Ising energy exhibits a critical behavior compatible 
with the sublattice order parameter at the low temperature region while the 
broad peaks in the sublattice susceptibilities and the lattice specific heat for 
0.05Hk are separate into the three peaks at the simulations with the O.OIHk 
heating rate. Consequently, there are a peak corresponding to each phase 
transitions for the successive AB (type-I)— F — AB (type-II)— F phase tran- 
sitions. On the other hand, the calculations for different lattice sizes ( L = 6, 
9 and 15) show that the model exhibits the AB (type-II) phase which occurs 
independently from linear dimension of lattice. 

At phase boundary between the AB (type-I) and the AB3(f) ground state 
regions, the model exhibits the multi AB (type-1)— F — AB (type-1)— F — P 
phase transition. In Fig. 8, the temperature variations of the thermodynamic 
quantities are shown for d — —15.8 parameter value. The ground state 
ordering for this parameter value is AB (type-I). At low temperature region, 
the sublattice order parameters TJia = mi, = Qa = Qb decreases sharply from 
unity to a value (~ 0.6) of ferromagnetic ordering and coincide with the 

— ma — Qc — Qd parameters at a transition temperature kTa/ J. For the 
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first order AB (type-I)— F phase transition, the sublattice order parameters 
and the lattice Ising energy are in s-shape nature (Fig. 8(a) and 8(c)) while 
the sublattice susceptibility and the lattice specific heat show a jump at the 
transition temperature (Fig. 8(b) and 8(d)). At a certain temperature above 
kTti/J , the sublattice order parameters part equally from nia — Qa — 0.6. 
In this case, the model has the AB (typc-1) ordering above a transition 
temperature kTf2/ J ■ At this transition, sublattice order parameters and the 
lattice Ising energy appear discontinuously and sublattice susceptibilities and 
the lattice specific heat exhibit a jump at the transition temperature kTt2/ J- 
As the temperature increases, the sublattice order parameters coincide each 
other at kTci/J and this caused the ferromagnetic ordering. At a high 
temperature above kTci/J, the sublattice order parameters nia decrease to 
the value {nia — 0) oi P ordering near the critical temperature kTc2/J- In 
the F — P phase transition, the sublattice susceptibilities and the lattice 
specific heat has a characteristic peak at kTc2/ J- 

In the interval —24 < d < —16, the model has the ABs{f) ground state 
ordering and shows the successive ABs{f) — F — P phase transitions through 
the k = —3 line. For a selected d = —20 value, the temperature variation 
of the thermodynamic quantities are illustrated in Fig. 9. In the AB^lf) 
ordering structure, three of the sublattices are occupied by B {Si — +1 
or —1) species { nib — rric — ''^a — ^, Qb — Qc — Qd — 1); while the other 
sublattice is occupied by A {Si = 0) species {nia = , = 0). As the 
temperature increases, — nic — and Qb — Qc — Qd decreases sharply 
from unity to a value (~ 0.7) of the F ordering and coincide with nia and Qa 
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at a transition temperature kTt/ J. However, the model exhibit the second 
order F — P phase transition at a critical temperature kTc/J above the 
kTt/ J. For the d = —20 value AB^i^f) — F and the F — P phase transitions 
are of the first and the second order, respectively. Except the d — —20 
value, the AB3{f) — F phase transition is of the second order in the interval 
-24 < d < -16. 

2.2. The phase diagrams 

The critical behavior of the BEG model near the phase boundaries is in- 
vestigated using CA algorithm through the k — —3 and d — —A lines. For 
k = —3, the estimated {kTc/J, d) phase diagram is shown in Fig. 10. As 
it is seen in Fig. 1, in the interval —24 < d < 0, the model demonstrates 
ferromagnetic (F), quadrupolar (Q) and staggered quadrupolar (SQ) order- 
ing regions compatible with the ground state phase diagram on the {k, d) 
plane through the k = —3 line. For k = —3, there are the different types of 
successive phase transitions as ABs{f) — F — P , the AB (type-I) —F — AB 
(type-II)-F - P, the AB (type-I) -F - P, the A3B(a)-A3B(f) - F-P 
and the AsB{a)—P. For all parameter values, the F — P phase transitions 
are of the second order. Thus, the paramagnetic P and the ferromagnetic 
F regions are separated with the second order critical line. This behavior 
is in compatible with the results of the other calculations [28 — 36]. On the 
other hand, the phase transitions from SQ orderings to F and P orderings 
are usually of the second order except for a few parameter values. The SQ 
ordering and other orderings are separated with second order critical line ac- 
cording to MFA results and CVM results [14, 21, 26], but this critical line is 
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of the first order according to CVM results [24] for k — —3. The CA results 
show that the line between the SQ and other orderings is not completely a 
second order critical line, because the model exhibits first order phase transi- 
tions for a few parameter value at the critical line. Consequently, our results 
roughly agree with the MFA ones for type of the critical lines. However, 
the previous works indicated that the SQ region generally separated a few 
of the sub-ordering regions as the ferrimagnetic, the antiquadrupolar and 
the dense ferromagnetic [14, 24, 26]. For the fee lattice, there are the Cu-Au 
iAsB{a), A^B{f), AB{f) and AB3{f) ) type orderings [24] in the SQ region 
of the phase diagram. However, the Cu-Au type systems exhibit also the 
AB^type-lT) ordering according to the experimental results [39, 41 — 47] and 
the theoretical studies [31,36,37]. Our calculations demonstrate that the 
model has the AB (type-II) ordering regions in addition to AsB{a), AsB(f), 
AB ( /) and AB^{f) ordering regions. It is seen from the {kTc/J, d ) phase 
diagram that there is the A^B{a) ordering in the interval — 8 < c? < 0, the 
74S(type-I) ordering in the interval —16 < d < — 8 and the AB^{f) ordering 
in the interval —24 < d < —16, at kT/ J — 0. However, in high temperatures, 
there arise the A2B[f) ordering region between the y4i?(type-l) and AB^{f) 
ordering regions and ^45 (type-II) ordering regions over the AB {type-1) or- 
dering region. At the phase boundary, the successive AB (type-I)— F — AB 
(type-I)— F — AB (type-II)— F — P phase transitions occur between the AB 
(type-I) and the AB:-i{f ) ground state regions at li = —15.9. The low tem- 
perature phase transitions at this parameter are of the first order. As a 
result of these phase transitions, two ground state orderings separate with 
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the first order phase transition hne. The phase diagram near this parameter 
has a cavity including ferromagnetic F ordering. As a result, the SQ region 
splits A^B{a ), AsB{f), A5(type-I), AB{type-II) and AB^^if) orderings in 
the different intervals of the d and the temperature values. The critical tem- 
peratures of AB (type-I)— F — AB (type-II) phase transitions shown as a 
band (thick line) because of being close each other on the phase diagram. 

The {kTc / J, k) phase diagram through the d = —A line is presented in 
Fig. 11. Through the d — —A line, our calculations indicate that, there 
are the ABs{f) ordering in the interval —15 < k < —1.35, the AB (type-I) 
ordering in the interval —2 < k < —1.5 and the ^3^(0) ordering for k < —2 
at kT/J — 0. The {kTc/ J, k) phase diagram for c? = — 4 does not have 
the AsB{f) ordering which occur between AS (type-I) and AsB{a) regions 
through the A; = — 3 hne. It indicates that the AsB{f) ordering region 
becomes narrow with increasing k value on the SQ ordering region of the 
{k, d) plane. In the interval —2<d< —1.5 above the AB{ty^Q-l) ordering, 
there occur AS (type-II) ordering regions for the certain parameter regions 
as the {kTc/ J, d) phase diagram for k — —3. The AB (type-II) ordering 
vanishes near the A'^B{P) ground state ordering region in the interval —2 < 
k < —1.75. For the different k values through the d — —4 line, the model 
exhibits the AB (type-I) -F- AS (type-II) -F-P, the AB (type-I) -F-P, 
the AsB{a)—F — P and the AS3(/) — F type successive phase transitions. 
In this phase transition, the A3B{a)—F, the A-iB{a)—P, the AB (type-1)— F 
phase transition at k — —1.8 value and the AB^{f)—F phase transitions in 
the —1.41 < k < —1.35 parameter region are of the first order while the 
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other transitions are of the second order. Consequently, the CA results for 
the {kTc/J, d) and the {kTc/J, k) phase diagrams through the k = —3 and 
d = —4 lines show that the model has the A^B, AB and ground state 
ordering regions on the (/c, 0?) plane which is illustrated in Fig. 1. However, 
there are different type critical behaviors in the phase boundaries between 
these regions and new ordering structures in the high temperature region 
over the ground state. 
3. Conclusion 

The spin-1 Ising BEG model has been simulated using the heating algo- 
rithm of the cellular automaton for the face-centered cubic lattice. The simu- 
lations show that the BEG model has the Cu-Au type stoichiometric A^B{a), 
AsB{f), AB{type-T), AB{type-ll) and ABs{f) orderings in the staggered 
quadrupolar region for k < —1. Furthermore, the AsB{a), AsB{f), AB(type- 
I) and y4i?3(/) ordering regions exhibit the different successive phase tran- 
sitions. The phase diagrams on the {kTc/J , d) and the {kTc/J, k) planes 
are obtained for k — —3 and d — —4, respectively. The cellular automaton 
(CA) results reconstruct the phase boundaries on the {kTc/J, d) phase di- 
agrams which indicated by the other calculations [14,24,26]. It is seen that 
the model exhibits the successive and the multi phase transitions from the 
first or the second order near the phase boundaries. The SQ ordering and the 
other orderings are separated with the second order critical line according to 
MFA results [14], but this critical line is of the first order according to CVM 
results [24] for k — —3. The CA results roughly agree with the MFA ones 
for type of the critical hues. In addition, it is seen that the fee BEG model 
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has the modulated AB (type-II) phase as indicated by the experimental re- 
sults [39, 41 — 47]. The modulated AB (typc-II) phase is obtained above the 
AB (type-I) phase which was indicated by the CVM [24] on the {kTc/J, 
d) plane. Furthermore, the heating rate is effective for separating the phase 
transitions which are close. The obtained results show that the BEG model 
is very suitable for studying with the binary alloys. Furthermore, the heating 
algorithm improved from Creutz cellular automaton is successful for exposing 
the critical behavior of the BEG model. 
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Figure Captions 

Fig. 1. The ground state phase diagram {k, d) of the BEG model on the 
fee lattice. 

Fig. 2. The four sublattices for an fee lattice. The circle, the square, 
the diamond and the triangle indicate the ma, m;,, rUc and rud sublattices, 
respectively. 

Fig. 3. For d — —4, the temperature dependence of (a) the lattice 
order parameters (M, Q) and sublatticc order parameters {ma,qa), (b) the 
sublattice susceptibility {xj, (c) the lattice Ising energy {Hj), and (d) the 
lattice specific heat (C/k) a,t k — —3 on L— 9. 

Fig. 4. For d — —7, the temperature dependence of (a) the sublattice 
order parameters {ma,qa), (b) the lattice susceptibility (x) and the sublat- 
tice susceptibility (x^), (c) the lattice Ising energy {Hj), and (d) the lattice 
specific heat (C/k) a,t k — —3 on L— 9. 

Fig. 5. For d — —7.7, the temperature dependence of (a) the sublattice 
order parameters {ma,qa), (b) the sublattice susceptibility {xj: (c) the lat- 
tice Ising energy (Hj) and (d) the lattice specific heat (C/k) a,t k — —3 on 
L= 9. 
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Fig. 6. For d — —12, the temperature dependence of (a) the sublattice 
order parameters {ma,qa), (b) the sublattice susceptibihty ixji (c) the lat- 
tice Ising energy (if/) and (d) the lattice specific heat {C /k) aX k = —3 on 
L= 9. 

Fig. 7. For d — —14.5, the temperature dependence of (a) the sublat- 
tice order parameters {rricqa), (b) the sublattice susceptibihty (c) the 
lattice Ising energy {Hj) and (d) the lattice specific heat iC/k) on L= 9 at 
k = -3. 

Fig. 8. For d — —15.8, the temperature dependence of (a) the sublattice 
order parameters {ma,qa), (b) the sublattice susceptibility {xj, (c) the lat- 
tice Ising energy {Hj) and (d) the lattice specific heat {C/k) at k — —3 on 
L= 9. 

Fig. 9. For d — —20, the temperature dependence of (a) the sublattice 
order parameters {ma,qa), (b) the sublattice susceptibility (c) the lat- 
tice Ising energy {Hj) and (d) the lattice specific heat {C/k) at k = —3 on 
L= 9. 

Fig. 10. The phase diagram in the {kTc/zJ, d) plane for /c = —3 on 
L= 9. The closed and open symbols indicate the second and the first order 
phase transitions, respectively. The thick line indicates the ferromagnetic ( 
F) ordering. 

Fig. 11. The phase diagram in the {kTc/zJ, k) plane for d — —A on 
L= 9. The closed and open symbols indicate the second and the first order 
phase transitions, respectively. The thick line indicates the ferromagnetic ( 
F) ordering. 
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